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Abstract 

The Supersymmetry Les Houches Accord (SLHA) provides a universal 
set of conventions for conveying spectral and decay information for su- 
persymmetry analysis problems in high energy physics. Here, we propose 
extensions of the conventions of the first SLHA to include various gener- 
alisations: the minimal supersymmetric standard model with violation of 
CP, R-parity, and flavour, as well as the simplest next-to-minimal model. 
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1 Introduction 



Supersymmetric (SUSY) extensions of the Standard Model rank among the most promising 
and well-explored scenarios for New Physics at the TeV scale. Given the long history of 
supersymmetry and the number of people working in the field, several different conventions 
for defining supersymmetric theories have been proposed over the years, many of which have 
come into widespread use. At present, therefore, no unique set of conventions prevails. In 
principle, this is not a problem. As long as everything is clearly and consistently defined, 
a translation can always be made between two sets of conventions. 

However, the proliferation of conventions does have some disadvantages. Results ob- 
tained by different authors or computer codes are not always directly comparable. Hence, 
if author/code A wishes to use the results of author/code B in a calculation, a consistency 
check of all the relevant conventions and any necessary translations must first be made - a 
tedious and error-prone task. 

To deal with this problem, and to create a more transparent situation for non-experts, 
the original SUSY Les Houches Accord (SLHAl) was proposed [1]. This accord uniquely 
defines a set of conventions for supersymmetric models together with a common interface 
between codes. The most essential fact is not what the conventions are in detail (they largely 
resemble those of [2]), but that they are consistent and unambiguous, hence reducing the 
problem of translating between conventions to a linear, rather than a factorial, dependence 
on the number of codes involved. At present, these codes can be categorised roughly as 
follows (see [3,4] for a review and on-line repository): 

• Spectrum calculators [5-8], which calculate the supersymmetric mass and coupling 
spectrum, assuming some (given or derived) SUSY-breaking terms and a matching to 
known data on the Standard Model parameters. 

• Observables calculators [9-19]; packages which calculate one or more of the fol- 
lowing: collider production cross sections (cross section calculators), decay partial 
widths (decay packages), relic dark matter density (dark matter packages), and indi- 
rect/precision observables, such as rare decay branching ratios or Higgs/electroweak 
observables (constraint packages). 

• Monte-Carlo event generators [20-28], which calculate cross sections through explicit 
statistical simulation of high-energy particle collisions. By including resonance decays, 
parton showering, hadronisation, and underlying-event effects, fully exclusive final 
states can be studied, and, for instance, detector simulations interfaced. 

• SUSY and CKM fitting programs [29-32] which fit model parameters to collider-type 
data. 

At the time of writing, the SLHAl has already, to a large extent, obliterated the need for 
separately coded (and maintained and debugged) interfaces between many of these codes. 
Moreover, it has provided users with input and output in a common format, which is more 
readily comparable and transferable. Finally, the SLHA convention choices are also being 
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adapted for other tasks, such as the SPA project [33]. We beheve, therefore, that the 
SLHA project has been useful, solving a problem that, for experts, is trivial but frequently 
occurring and tedious to deal with, and which, for non-experts, is an unnecessary headache. 

However, SLHAl was designed exclusively with the MSSM with real parameters and 
R-parity conservation in mind. Some recent public codes [6,7,17,18,34-38] are either im- 
plementing extensions to this base model or are anticipating such extensions. It therefore 
seems prudent at this time to consider how to extend SLHAl to deal with more general 
super symmetric theories. In particular, we will consider R-parity violation (RPV), flavour 
violation (FLV), and CP- violating (CPV) phases in the minimal supersymmetric standard 
model (MSSM). We will also consider next-to- minimal models (i.e., models in which the 
MSSM field content is augmented by a gauge-singlet chiral superfield) which we shall col- 
lectively label by the acronym NMSSM. 

Rather than giving exhaustive historical references for all concepts used in this article, 
we provide a list of useful and pedagogical reviews to whose contents and references in turn 
we refer. For the various topics treated in the article, these reviews are: 

• SUSY [39], 

• FLV [40], 

• Neutrinos [41], 

• RPV [42,43], 

• CPV [44], 

• NMSSM [45], 

• SUSY Tools [3,4]. 

There is clearly some tension between the desirable goals of generality of the models, 
ease of implementation in programs, and practicality for users. A completely general accord 
would be useless in practice if it was so complicated that no one would implement it. We 
have agreed on the following for SLHA2: for the MSSM, we will here restrict our attention 
to either CPV or RPV, but not both. For RPV and flavour violation, we shall work in the 
Super-CKM/PMNS basis, as defined in sections O and H^l For the NMSSM, we define 
one catch-all model and extend the SLHAl mixing only to include the new states, with CP, 
R-parity, and flavour still assumed conserved. 

To make the interface independent of programming languages, compilers, platforms etc, 
the SLHAl is based on the transfer of three different ASCII files (or potentially a character 
string containing identical ASCII information): one for model input, one for spectrum 
calculator output, and one for decay calculator output. We believe that the advantage of 
implementation independence outweighs the disadvantage of codes using SLHAl having to 
parse input. Indeed, there are tools to assist with this task [46-48]. 

Care was taken in SLHAl to provide a framework for the MSSM that could easily be 
extended to the cases listed above. The conventions and switches described here are de- 
signed to be a superset of those of the original SLHAl and so, unless explicitly mentioned in 
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the text, we will assume the conventions of the original accord [1] implicitly. For instance, 
all dimensionful parameters quoted in the present paper are assumed to be in the appro- 
priate power of GeV, all angles are in radians, and the output formats for SLHA2 data 
BLOCKS follow those of SLHAl. In a few cases it will be necessary to replace the original 
conventions. This is clearly remarked upon in all places where it occurs, and the SLHA2 
conventions then supersede the SLHAl ones. 

2 Extensions of SLHAl 

Since its first publication, a few useful extensions to the SLHAl have been identified. These 
are collected here for reference and are independent of the more general SUSY models 
discussed in subsequent sections. (Also note the recent proposal for a joint SLHA+LHEF 
format for BSM event generation [49,50].) 

Firstly, we introduce additional optional entries in the SLHAl block EXTPAR to allow 
for using either the A'^ or pole masses as input instead of the parameter m^(Minput) 
defined in [1]. 

Secondly, to allow for different parameters to be defined at different scales (e.g., n defined 
at Mewsb, the remaining parameters defined at Mjnput) we introduce a new optional block 
QEXTPAR which, if present, overrides the default MINPAR and EXTPAR scale choices for specific 
parameters, as defined below. 

While there is no obligation on codes to implement these extensions, we perceive it as 
useful that the accord allows for them, enabling a wider range of input parameter sets to be 
considered. The entries defined in EXTPAR and QEXTPAR in the SLHA2 are thus (repeating 
unchanged EXTPAR entries for completeness): 

BLOCK EXTPAR 

Optional input parameters for no n-minimal/non- universal models. This block may be en- 
tirely absent from the input file, in which case a minimal type of the selected SUSY breaking 
model will be used. When block EXTPAR is present, the starting point is still a minimal 
model with parameters as given in MINPAR [1] but with each value present in EXTPAR re- 
placing the minimal model value of that parameter, as applicable. If MINPAR is not present, 
then all model parameters must be specified explicitly using EXTPAR. All scale-dependent 
parameters are understood to be given in the DR scheme. 
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Input scale 



: Minpuf Input scale for EXTPAR entries in SUGRA, AMSB, and general 
MSSM models. If absent, the GUT scale derived from gauge unification 
will be used as input scale. Note that this parameter has no effect in 
GMSB scenarios where the input scale by definition is identical to the 
messenger scale, Mmess- A special case is when Q = Mewsb = y/iWj^ 
is desired as input scale, since this scale is not known beforehand. This 
choice can be invoked by giving the special value Mjnput = — 1- To 
define an alternative input scale for one or more specific parameters, see 
QEXTPAR below. 



Gaugino Masses 



1 

2 
3 



Mi(Minput)- U{1)y gaugino (Bino) mass. 
M2(Minput). SU{2)l gaugino (Wino) mass. 
M3(Minput). SU{?,)c gaugino (gluino) mass. 



Trilinear Couplings 



11 


: A(Minput) 


Top trilinear coupling. 


12 


: Afe(Minput) 


Bottom trilinear couplin 


13 


: ^r(Mi„p,t) 


. Tau trilinear coupling. 



Higgs Parameters 

— Only one of the parameter sets mj^J, {fi,m\), (/i, mAo), or (/i, mH+) should be 

given, they merely represent different ways of specifying the same parameters. 

21 : (Minput). Down type Higgs mass squared. 

22 : m^^^lMnput)- Up type Higgs mass squared. 

23 : /x(Minput). /i parameter. 

24 : m^(Minput)- Tree-level pseudoscalar Higgs mass parameter squared, as 

defined in [1]. 
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25 : tan/3(Minput)- If present, this value of tan/3 overrides the one in MINPAR, 

and the input scale is taken as Mi^put rather than m^. 

26 : m^o. Pseudoscalar Higgs pole mass. May be given instead of m^(Minput) • 

27 : mH+. Charged Higgs pole mass. May be given instead of m^(Minput). 



Sfermion Masses 


31 


■ '"'e^V input 


Tipft 1 ^^P'pn ^fril^ir Ipnton ma^^ 

' M , 1 1 1 _L CLV^AJ.. iJ\jCAiL(JLl. AV^L/ljWAJ. AXACAjOO« 


32 


■ m- (M- A 
• Mi V input J 


Tipft y^^^^ffpn ^f/^li^r Ipnton ma^^ 

' 1 1 1 1 ^ CLV^IJ.. 0*-jCAj1C(iX 1V>1_/UWJ.J. AXACAjOO" 


33 


' Tn - ( A/f;„„.,^ 1 
• \^ J input J ' 


Tjpft sralar lenton mass 




■ '"■e^ I, -'"input J • 






■ ''i^ifl I^J^^mput J 


R 1 cr n"!' cr*?! 1 T* Tnn on Tn 

± Lltiliu ol^CLlCLl iilLLvJii iiiCLoo. 


36 


• '^fij(^input)- 


Right scalar tau mass. 


41 


• '^giL ("^input) 


. Left I'^^gen. scalar quark mass. 


42 


• '^g2L ("^input) 


. Left 2"'^gen. scalar quark mass. 


43 


• '^gsL ("^input) 


. Left 3''*^gen. scalar quark mass. 


44 


• '^Mfl(^input) 


Right scalar up mass. 


45 


: mc^(Minput)- 


Right scalar charm mass. 


46 


: mj^(Minput). 


Right scalar top mass. 


47 


: mj^(Minput)- 


Right scalar down mass. 


48 


: m^-^(Minput)- 


Right scalar strange mass. 


49 


: m^^(Minput)- 


Right scalar bottom mass. 



Other Extensions 

51 : Ni (GMSB only). f/(l)y messenger index (defined as in ref. [51]). 

52 : N2 (GMSB only). SU{2)l messenger index (defined as in ref. [51]). 

53 : iVs (GMSB only). SU{'i)c messenger index (defined as in ref. [51]). 
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BLOCK QEXTPAR 

Optional alternative input scales for specific parameters. This block should normally be 
absent, in which case the default input scale or Minput (see EXTPAR 0) will be used for all 
parameters. We stress that most codes cannot be expected to allow for multiple arbitrary 
scale choices, so the relevant manual and output should be carefully checked to make sure 
the desired behaviour is obtained. Currently defined entries are: 



-1 
i 


Qmi- 


Input scale for Mi. 


o 
A 




Input scale for M2. 


o 
O 


Qms- 


Input scale for M3. 


1 1 


QAu- 


Input scale for up-type squark trilinear couplings. 


12 




Input scale for down-type squark trilinear couplings. 


13 


Qac 


Input scale for charged slepton trilinear couplings. 


21 




. Input scale for m^^ . 


99 




. Input scale for m'jj^. 


23 


Q^. Input scale for /i. 


24 




Input scale for m^, as defined in [1]. 


25 


Q tan 13 


. Input scale for tan (3. 


31 




Input scale for all left-handed slepton mass terms. 


34 




. Input scale for all right-handed slepton mass terms. 


41 




. Input scale for all left-handed squark mass terms. 


44 




. Input scale for all right-handed up-type squark mass terms. 


47 


Qrrij 


. Input scale for all right-handed down-type squark mass terms. 



3 Model Selection 

To define the general properties of the model, we propose to introduce global switches in 
the SLHAl model definition block MODSEL, as follows. Note that the switches defined here 
are in addition to the ones in [1]. 

BLOCK MODSEL 

Switches and options for model selection. The entries in this block should consist of an 
index, identifying the particular switch in the listing below, followed by another integer or 
real number, specifying the option or value chosen: 
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(Default=0) Choice of particle content. Switches defined are: 

: MSSM. This corresponds to SLHAl. 

1 : NMSSM. The blocks defined in section [5] should be present. 



(Default=0) R-parity violation. Switches defined are: 

: R-parity conserved. This corresponds to the SLHAl. 

1 : R-parity violated. The blocks defined in section should be 

present. 

(Default=0) CP violation. Switches defined are: 

: CP is conserved. No information even on the CKM phase is 

used. This corresponds to the SLHAl. 

1 : CP is violated, but only by the standard CKM phase. All 

other phases are assumed zero. 

2 : CP is violated. Completely general CP phases allowed. Imag- 

inary parts corresponding to the entries in the SLHAl block 
EXTPAR can be given in IMEXTPAR (together with the CKM 
phase). In the case of additional SUSY fiavour violation, 
imaginary parts of the blocks defined in section 14.11 should 
be given, again with the prefix IM, which supersede the cor- 
responding entries in IMEXTPAR. 

(Default=0) Flavour violation. Switches defined are: 

: No (SUSY) fiavour violation. This corresponds to the SLHAl. 

1 : Quark fiavour is violated. The blocks defined in section 14.11 

should be present. 

2 : Lepton fiavour is violated. The blocks defined in section 14.11 

should be present. 

3 : Lepton and quark fiavour is violated. The blocks defined in 

section H]T] should be present. 



4 General MSSM 

For convenience, we here repeat the definitions of the field content, superpotential and soft 
SUSY-breaking potential of the MSSM in the notation of [1]. 

Specifically, the chiral superfields of the MSSM have the following SU{3)c ® SU{2)l ® 
U{1)y quantum numbers 



L: (1,2,-i), E: (1,1,1), g:(3,2,i), t/:(3,l,-|), 

D : (3,l,i), iJi: (1,2,-i), (i,2,i) , (1) 

the superpotential (omitting RPV terms, see section 14.21) is written as 

Wmssm = eaA{YE)^,H'lL\E^ + {YD)i,HlQ\Dj + {Yu)^,Hl^^^^ , (2) 
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and the trilinear and bilinear soft SUSY-breaking potentials V3 and V2 are 

Vs = ea, E [iTEhHfLl~e]^ + {TnhmQ\J]^ + {Tu).,HIQIu*^] + h.c. , (3) 

Ui^ijnDijU*^ + + ^iRijnldij^]n - {mleabH^H^ + h.c.) , (4) 

where a tilde over the symbol for a quark or lepton superfield denotes its scalar component 
(note however that we define, e.g., mJj as the scalar component of U). Throughout this 
section, we denote SU{2)l fundamental representation indices by a, 6 = 1,2 and genera- 
tion indices by i,j = 1,2,3. Colour indices are everywhere suppressed, since only trivial 
contractions are involved, eab is the totally antisymmetric tensor, with ei2 = e^^ = 1. 

4.1 Flavour Violation 

4.1.1 The quark sector and the super-CKM basis 

Within the MSSM there are in general new sources of flavour violation arising from a 
possible misalignment of quarks and squarks in flavour space. The severe experimental 
constraints on flavour violation have no direct explanation in the structure of the uncon- 
strained MSSM which leads to the well-known supersymmetric flavour problem. 

The Super-CKM basis of the squarks is very useful in this context because in that basis 
only physically measurable parameters are present. In the Super-CKM basis the quark mass 
matrix is diagonal and the squarks are rotated in parallel to their superpartners. Actually, 
once the electroweak symmetry is broken, a rotation in flavour space 

D° = VdD, U° = VuU, f/*5, U°= U:U, (5) 

of all matter superfields in the (s) quark superpotential 

Wq = eat [(Id),, H^^Q',"DJ + (Yu),^ H'.QrUf ] (6) 

brings fermions from the interaction eigenstate basis {d1,u1,d°j^,u°p,} to their mass eigen- 
state basis {dL,UL,dR,UR}: 

dl = Vddi , = VuUl , d°j^ = Uddn , = UuUr , (7) 

and the scalar superpartners to the basis {^l, n^, c//?, £t/j}. Through this rotation, the 
Yukawa matrices and Yjj are reduced to their diagonal form Yd and Yu: 

(Yoh = (U^YEVdh = V2^ , (Yu)u = iUlY^^V^)^ = . (8) 

Tree-level mixing terms among quarks of different generations are due to the misalignment 
of Vd and Vu, expressed via the CKM matrix 

VcKM = V^Vd , (9) 
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which is proportional to the tree-level undLjW^ , uudjijH^ , and URidLjH'^ couplings (i, j = 
1, 2, 3). This is also true for the supersymmetric counterparts of these vertices, in the limit 
of unbroken supersymmetry. 

In the super-CKM basis the 6x6 mass matrices for the up-type and down-type squarks 
are defined as 

where = (ul, cl, ti, £t/j, c/j, t/j)^ and = {dL,SL,bL,dR,SR,bR)'^. We diagonalise the 
squark mass matrices via 6x6 unitary matrices Ru,d, such that Ru^d-M^^Rl d diagonal 
matrices with increasing mass squared values. The flavour-mixed mass matrices read: 

, VcKM ml ^c^KM + ml + Dull - fim^ cot B 

Mi = t „ , 11 



^Tu - yu*m„ cot p mi + mi + Du 



RR 



^ ( ^Q + '^d + DdLL ^U- fxnidtanP \ ^^^^ 

In the equations above we introduced the 3x3 matrices 

rnl^VlmlVd, ml^Ulmfu^, ih\^U\mf Ud. (13) 

f^^[/tr^K, %^\j\TlVd. (14) 

where the un-hatted mass matrices vp?q ^ ^ and trilinear interaction matrices Tu d are given 
in the interaction basis. 

The matrices ra^^d are the diagonal up-type and down- type quark masses and DfLL,RR 
are the D-terms given by: 

Df LL,RR = COS 2P m| {Tf-Qf sin^ Ow) h , (15) 

which are also flavour diagonal. Here, Qf is the electric charge of the left-handed chiral 
supermultiplet to which the squark belongs, i.e., it is 2/3 for U and —2/3 for W^. Note 
that the up-type and down-type squark mass matrices in eqs. (!TT|) and (fT2|) cannot be 
simultaneously flavour-diagonal unless ttiq is flavour- universal (i.e., proportional to the 
identity in flavour space). 



4.1.2 The lepton sector and the super-PMNS basis 

For the lepton sector, we adopt a super-PMNS basis, as defined in this section. 

Neutrino oscillation data have provided a strong indication that neutrinos have masses 
and that there are flavour-changing charged currents in the leptonic sector. One popular 
model to produce such effects is the see-saw mechanism, where right-handed neutrinos have 
both Majorana masses as well as Yukawa couplings with the left-handed leptons.When the 
heavy neutrinos are integrated out of the effective field theory, one is left with three light 
approximately left-handed neutrinos which are identified with the ones observed experi- 
mentally. There are other models of neutrino masses, for example involving SU(2) Higgs 
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triplets, that, once the triplets have been integrated out, also lead to effective Majorana 
masses for the neutrinos. Here, we cover all cases that lead to a low energy effective field 
theory with Majorana neutrino masses and one sneutrino per family. In terms of this low 
energy effective theory, the lepton mixing phenomenon is analogous to the quark mixing 
case and so we adapt the conventions defined above to the leptonic case. 

After electroweak symmetry breaking, the neutrino sector of the MSSM contains the 
Lagrangian pieces (in 2-component notation) 

C = -^u°^{m^)u° + h.c., (16) 

where m^, is a 3 x 3 symmetric matrix. The interaction eigenstate basis neutrino fields 
are related to the mass eigenstate ones u by 

= V,u, (17) 

reducing the mass matrix to its diagonal form rhy 

{rhy)ii = {VjmyVy)ii = m^^. (18) 

The charged lepton fields have a 3 x 3 Yukawa coupling matrix defined in the (s)lepton 
superpotential 

We = eab{YE),jHlL^°E°, (19) 

where the charged lepton interaction eigenstates {e^, e^} are related to the mass eigenstates 
{eL,eR, } by 

el = VeeL and e?^ = f/eC^j. (20) 
The equivalent diagonalised charged lepton Yukawa matrix is 

iYE)u={UlY^V;)u = V2^ . (21) 

Lepton mixing in the charged current interaction can then be characterised by the PMNS 
matrix 

UpMNS = V}V, , (22) 

which is proportional to the tree-level eLiVjW~ and emVjH' couplings (z, j = 1, 2, 3). This 
is also true for the supersymmetric counterparts of these vertices, in the limit of unbroken 
sup er sy mmet r y. 

Rotating the interaction eigenstates of the sleptons identically to their leptonic coun- 
terparts, we obtain the super-PMNS basis for the charged sleptons and the sneutrinos, 
described by the Lagrangiar0 

= -^iMl^e - ^iMl^,, (23) 

^We here neglect the possible term i>^A^~$j/. 
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where = {ly^, u^, u^-Y and $e = (cl, /ii, t/,, e/j, /ir, t/?)^. M.\ is the 6x6 matrix 

/ m| + + -^f^ - /ime tan/5 \ 

7W- = ^ , (24) 

V ff^E - yu*me tan/5 m| + + D^^^ y 

and A^l is the 3x3 matrix 

-^^ = UImns rril Upmns + D^^l. (25) 

where -DgLL and are given in eq. (fT5]) . In the equations above we introduced the 3x3 
matrices 

ml = V} ml Ve, ml = Ul mf Ue , (26) 

fE = UlT^Ve, (27) 

where the un-hatted mass matrices m| ^ and the trihnear interaction matrix Te are given 
in the interaction basis. We diagonahse the charged slepton and sneutrino mass matrices 
via the unitary 6x6 and 3x3 matrices Re^u respectively. Thus, Re^u-Mg^uRl^u diagonal 
with increasing entries toward the bottom right of each matrix. 



4.1.3 Explicit proposal for SLHA2 

As in the SLHAl [1], for all running parameters in the output of the spectrum file, we 
propose to use definitions in the modified dimensional reduction (DR) scheme. The basis 
is the super-CKM/PMNS basis as defined above, that is the one in which the Yukawa 
couplings of the SM fermions, given in the DR scheme, are diagonal. Note that the masses 
and vacuum expectation values (VEVs) in eqs. ([H]), (flSll . and (1211) must thus also be the 
running ones in the DR scheme. 

The input for an explicit implementation in a spectrum calculator consists of the fol- 
lowing information: 

• By default, all input SUSY parameters are given at the scale Minput as defined in 
the SLHAl block EXTPAR (see above). In principle, advanced codes may also allow 
for separate input scales for the sfermion mass matrices and trilinear couplings, via 
the block QEXTPAR defined above, but we emphasise that this should be regarded as 
non-standard. 

• For the SM input parameters, we take the Particle Data Group (PDG) definition: 
lepton masses are all on-shell. The light quark masses m„ ^ are given at 2 GeV, and 
the heavy quark masses are given as mc{mc)^^ , mb(mb)^^ and m°°~^^'^^\ The latter 
two quantities are already in the SLHAl. The others are added to SMINPUTS in the 
following manner (repeating the SLHAl parameters for convenience): 

1 : a~^{mz)^^- Inverse electromagnetic coupling at the Z pole in the MS 

scheme (with 5 active fiavours). 

2 : Gp- Fermi constant (in units of GeV~^). 



13 



3 : as(m^)^^. Strong coupling at the Z pole in the MS scheme (with 5 





7\c\\\T(^ fli^AronTc; 1 


4 


: rriz, pole mass. 


5 


: mb(mh)'^^. h quark running mass in the MS scheme. 


6 


: m+, pole mass. 


7 


: m,-, pole mass. 


8 


: rriy^, pole mass. 


11 


: me, pole mass. 


12 


: rriy^, pole mass. 


13 


: m^, pole mass. 


14 


: m,^2, pole mass. 


21 


: mrf(2 GeV)'^^. c? quark running mass in the MS scheme. 


22 


: mM(2 GeV)^^. u quark running mass in the MS scheme. 


23 


: m<i(2 GeV)'^^. s quark running mass in the MS scheme. 


24 


: mc{mcy^^- c quark running mass in the MS scheme. 



The FORTRAN format is the same as that of SMINPUTS in SLHAl [1]. 

• VcKM- the input CKM matrix in the Wolfenstein parameterisatior^, in the block 
VCKMIN. Note that present CKM studies do not precisely define a renormalisation 
scheme for this matrix since the electroweak effects that renormalise it are highly 
suppressed and generally neglected. We therefore assume that the CKM elements 
given by PDG (or by UTFiT [31] and CKMFlTTER [32], the main collaborations 
that extract the CKM parameters) refer to SM MS quantities defined at Q = mz, to 
avoid any possible ambiguity. VCKMIN should have the following entries 

1 : A 

2 : A 

3 : p 

4 : f] 

The FORTRAN format is the same as that of SMINPUTS above. 

• ^PMNs: the input PMNS matrix, in the block UPMNSIN. It should have the PDG 
parameterisation in terms of rotation angles [52] (all in radians): 

1 : 6i2 (the solar angle) 

2 : $23 (the atmospheric mixing angle) 

^For the Wolfenstein parameters we use the PDG definition, eq.(11.4) of [52], which is exact to aU orders 
in A. 
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4 



3 



6 



5 



^13 (currently only has an upper bound) 
5i3 (the Dirac CP-violating phase) 
«! (the first Majorana CP-violating phase) 
(the second CP-violating Majorana phase) 



The FORTRAN format is the same as that of SMINPUTS above. Majorana phases 
have no effect on neutrino oscillations. However, they have physical consequences in 
the case of, for example, /3/90i/ decay of nuclei. 

• (rn|)°^, (m|)°-^, (m|)°^, (m|)°-^, (m|)°-^: the squark and slepton soft SUSY- 
breaking masses at the input scale in the super-CKM/PMNS basis, as defined above. 
They will be given in the new blocks MSQ2IN, MSU2IN, MSD2IN, MSL2IN, MSE2IN, with 
the FORTRAN format 

(lx,I2,lx,I2,3x,lP,E16.8,0P,3x, '#\lx,A) . 

where the first two integers in the format correspond to i and j and the double 
precision number to the soft mass squared. Only the "upper triangle" of these matrices 
should be given. If diagonal entries are present, these supersede the parameters in 
the SLHAl block EXTPAR. 

• (tV)?-^, (Td)?-^, and (Tg;)^^-^: the squark and slepton soft SUSY-breaking trilinear 
couplings at the input scale in the super-CKM/PMNS basis. They will be given in the 
new blocks TUIN, TDIN, TEIN, in the same format as the soft mass matrices above. If 
diagonal entries are present these supersede the A parameters specified in the SLHAl 
block EXTPAR [1]. 

For the output, the pole masses are given in block MASS as in SLHAl (note, however, that 
some PDG numbers have different assignments in SLHA2, see below) and the DR and 
mixing parameters as follows: 

• ('^Di^-^' ('^1)^-^' ('^i)!^-^' (^^1)°^: the squark and slepton soft SUSY- 
breaking masses at scale Q in the super-CKM/PMNS basis. Will be given in the 
new blocks MSQ2 Q=. . ., MSU2 Q=. . ., MSD2 Q=. . ., MSL2 Q=. . ., MSE2 Q=. . ., with 
formats as the corresponding input blocks MSX2IN above. 

• (T(/)i]-^! K^ofif^i and (Te)^]-^: The squark and slepton soft SUSY-breaking trilinear 
couplings in the super-CKM/PMNS basis. Given in the new blocks TU Q=. . ., TD 
Q=. . ., TE Q=. . ., which supersede the SLHAl blocks AD, AU, and AE, see [1]. 

• (n/)H^, ^dW. (^s)IF: the diagonal DR Yukawas in the super-CKM/PMNS basis, 
with Y defined by eqs. ([8]) and (pTj) . at the scale Q. Given in the SLHAl blocks YU 
Q=. . ., YD Q=. . ., YE Q=. . ., see [1]. Note that although the SLHAl blocks provide 
for off-diagonal elements, only the diagonal ones will be relevant here, due to the 
CKM/PMNS rotation. 
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The entries of the DR CKM matrix at the scale Q. The real and imaginary parts are 
given in VCKM Q= . . . and IMVCKM Q= . . . , respectively. The format of the individual 
entries is the same as for mixing matrices in the SLHAl. Note that the complete 
matrix should be output, i.e., all entries should be included. 



The entries of the DR PMNS matrix at the scale Q. The real and imaginary parts 
are given in UPMNS Q= . . . and IMUPMNS Q= . . . , respectively, with entries defined as 
for the VcKM output blocks above. 

The squark and slepton masses and mixing matrices should be defined as in the 
existing SLHAl, e.g. extending the t, h and f mixing matrices to the 6x6 case. More 
specifically, the new blocks Ru =USQMIX Rd =DSQMIX, R^ =SELMIX and the 3 x 3 
matrix for =SNUMIX specify the composition of the mass eigenstates in terms of 
the super-CKM/PMNS basis states according to the following definitions: 



/ 1000001 \ 
1000003 
1000005 
2000001 
2000003 

V 2000005 / 

/ 1000002 \ 
1000004 
1000006 
2000002 
2000004 

V 2000006 J 

( 1000011 \ 
1000013 
1000015 
2000011 
2000013 

V 2000015 ) 

1000012 
1000014 
1000016 



DSQMIXi 



mass— ordered 



USQMIXi 



mass— ordered 



SELMIX 



mass— ordered 



( ~dL \ 
Sl 
h 

dR 

V bR J 
( \ 

h 

Ur 
cr 

V tR J 

( eL \ 

tr 



(2^ 



super-CKM 



(29) 



super-CKM 



(30) 



\rR / 



super-PMNS 



SNUMIX 



(31) 



mass— ordered 



V 



super-PMNS 



Note! A potential for inconsistency arises if the masses and mixings are not calculated 
in the same way, e.g. if radiatively corrected masses are used with tree-level mixing 
matrices. In this case, it is possible that the radiative corrections to the masses shift 
the mass ordering relative to the tree-level. This is especially relevant when near- 
degenerate masses occur in the spectrum and/or when the radiative corrections are 
large. In these cases, explicit care must be taken especially by the program writing 
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the spectrum, but also by the one reading it, to properly arrange the rows in the order 
of the mass spectrum actually used. 



Optionally, we allow for the possibility of the scalar and pseudoscalar components 
of the sneutrinos to be treated separately. In this case, we define separate PDG codes 
and mixing matrices for the scalar and pseudoscalar sneutrinos, as follows: 



/ 1000012 \ 




1000014 




V 1000016 ) 




1 1000017 ^ 




1000018 


H 


V 1000019 ) 






SNSMIX,; 



mass— ordered 



SNAMIX 



mass— ordered 



V2Re{i>e} 

y2Re {vr} 
j V2Im{i>e} 
V V2Im{z>,} 



If present, SNSMIX and SNAMIX supersede SNUMIX. 



supcr-PMNS 



super-PMNS 



(32) 



(33) 



4.2 R-Parity Violation 

We write the R-parity violating superpotential in the interaction basis as 

W^RPV = ' 



+-K,ke.y^tD]Dl (34) 



■1 
2' 

1 

2' 

where x,y,z = 1, . . . , 3 are fundamental SU(3)c indices and e^yz is the totally antisymmetric 
tensor in 3 dimensions with eus = +1. In eq. flMl) . Xijk, X'ljk ki break lepton number, 
whereas A^'^-^, violate baryon number. To ensure proton stability, either lepton number 
conservation or baryon number conservation is usually still assumed, resulting in either 
\jk = Kjk = = or A-'^-fc = for all k = 1,2, 3. 

The trilinear R-parity violating terms in the soft SUSY-breaking potential are 

Vs.RPV = Cafe -{T)ijk^L^LKR+ {T')ijk^LQ]LdkR 

+\{T'%ke.y.urRdTRdlR + h.c. . (35) 

Note that we do not factor out the A couplings (e.g. as in Tijk/Xijk = Ax,ijk) in order to 
avoid potential problems with Xijk = but Tijk ^ 0. This usage is consistent with the 
convention for the R-conserving sector elsewhere in this report. 

The bilinear R-parity violating soft terms (all lepton number violating) are 

1^2,RPV = -ea.D.L'^^Hl + LUml^H^ + h.c. . (36) 

When lepton number is not conserved the sneutrinos may acquire vacuum expectation 
values (VEVs) {i'e,n,T) = Ve,fi,T/V2- The SLHAl defined the VEV v, which at tree level is 



equal to2mz/^Jg^ + g'^ ~ 246 GeV; this is now generalised to 

V = ^vf + v^ + v^ + vl + v^ . (37) 
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The addition of sneutrino VEVs allows for various different definitions of tan /3, but we here 
choose to keep the SLHAl definition tan/? = v^jvx. 

For input/output, we use the super-CKM/PMNS basis throughout, as defined in sec- 
tion 14.11 with the following considerations specific to the R-parity violating case. 

Firstly, the (i-quark mass matrices are given by 

V2{md)ij = {YD)ijVi + X'kijVk . (38) 

where Vk are the sneutrino VEVs. Secondly, in the lepton number violating case, the PMNS 
matrix can only be defined consistently by taking into account the 1-loop contributions 
induced by the lepton- number violating couplings (see, e.g., [43]). We here restrict our 
attention to scenarios in which there are no right-handed neutrinos and, thus, neutrino 
masses are generated solely by the lepton number violating couplings. In this case, the 
PMNS matrix is not an independent input but an output. 

For definiteness, and to keep the changes with respect to the R-parity conserving case as 
limited as possible, we define the super-CKM basis as the one where the Yukawa couplings 
Yd and Yu are diagonal. The PMNS basis is defined as the basis where Ye is diagonal and 
the loop-induced neutrino mass matrix is diagonalised. In this way one obtains a uniquely 
defined set of parameters: 

Ajjfc = \-st^p,riVe,sjU\^^j^ , (39) 
Kjk = Kst^y,TiVd,sjUltk , (40) 

ki = K^V^yj , (41) 

Kjk = \stUu,riUd,sjUd,tk ) (42) 

where the fermion mixing matrices are defined in section 14. 1[ The Lagrangian for the 
quark-slepton interactions then takes the following form: 

C = -X'ijk^idRkdLj + KskUpMNS,riVcKM,sjL,idRkULj + h.C. . (43) 

Similarly one obtains the soft SUSY breaking couplings in this basis by replacing the super- 
potential quantities in eqs. ( I39i) -( l42l) by the corresponding soft SUSY breaking couplings. 
In addition we define: 

= ^UrK^H, ■ (44) 



4.2.1 Input/Output Blocks 

As mentioned above, we use the super-CKM/PMNS basis throughout, for both superpoten- 
tial and soft SUSY-breaking terms. This applies to both input and outputE]. The naming 
convention for input blocks is BLOCK RV#IN, where the '#' character represents the name 

•^A code may need to convert internally the parameters to the interaction basis. In this case it must 
supply - or take as additional inputs - the individual rotation matrices of quark and lepton superfields 
entering eqs. (|5g)) - P^ . 
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of the relevant output block given below (thus, for example, the "LLE" couplings in the 
super-PMNS basis, Xijk, would be given in BLOCK RVLAMLLEIN). 

Default inputs for all R-parity violating couplings are zero. The inputs are given at scale 
Mnput, as described in SLHAl (again, if no Mnput is given, the GUT scale is assumed), 
and follow the output format given below (with the omission of Q= . . .). In addition, the 
known fermion masses should be given in SMINPUTS as defined in section I4.1.3[ 

The dimensionless super-CKM/PMNS couplings Xijk, Kjk^ ^"^^ Kjk given in BLOCK 
RVLAMLLE , RVLAMLQD , RVLAMUDD Q= . . . respectively. The output standard should cor- 
respond to the FORTRAN format 

(lx,I2,lx,I2,lx,I2,3x,lP,E16.8,0P,3x, ,lx,A) . 

where the first three integers in the format correspond toi, j, and k and the double precision 
number is the coupling. 

fijk, flj^, and fl'-^ are given in BLOCK RVTLLE, RVTLQD, RVTUDD Q= ... in the same 

format as for the A couplings above. 

The bilinear superpotential and soft SUSY-breaking terms ki, Di, and m| ^ and the 
sneutrino VEVs are given in BLOCK RVKAPPA, RVD, RVM2LH1, RVSNVEV Q= ! .] respec- 
tively, in the format 

(lx,I2,3x,lP,E16.8,0P,3x, ,lx,A) . 

The input and output blocks for R-parity violating couplings are summarised in Tab. [TJ 
As for the R-conserving MSSM, the bilinear terms (both SUSY-breaking and SUSY-respect- 
ing ones, including fi) and the VEVs are not independent parameters. They become related 
by the condition of electroweak symmetry breaking. Thus, in the SLHAl, one had the pos- 
sibility either to specify and rnjj^ or fi and m\. This carries over to the RPV case, 
where not all the parameters in the input blocks RV. . . IN in Tab. [1] can be given simulta- 
neously. Specifically, of the last 4 blocks only 3 are independent. One block is determined 
by minimising the Higgs-sneutrino potential. We do not here insist on a particular choice 
for which of RVKAPPAIN, RVDIN, RVSNVEVIN, and RVM2LH1IN to leave out, but leave it up to 
the spectrum calculators to accept one or more combinations. 

4.2.2 Particle Mixing 

In general, the neutrinos mix with the neutralinos. This requires a change in the definition 
of the 4x4 neutralino mixing matrix to a 7 x 7 matrix. The Lagrangian contains the 
(symmetric) neutrino/neutralino mass matrix as 



in the basis of 2-component spinors tp^ = {ue, Ufxji'T, —ib, —iw^ ,hi, h2)^ . We define the 
unitary 7x7 neutrino/neutralino mixing matrix A^ (block RVNMIX), such that: 




mass 



^V^°^A^^oV^° + h.c. , 



(45) 




(46) 



diag(?Ti-()) X' 
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Input block 


Output block 


data 


Vl\n AMT T TTTM 


vi\n AMT T tr 


i j k Xijk 


vi\n AMT nr>Ti\T 


vi\n AMT ni^ 




RVLAMUDDIN 


RVLAMUDD 


^ J k \jk 


RVTLLEIN 


RVTLLE 


i j k Tijk 


RVTLQDIN 


RVTLQD 




RVTUDDIN 


RVTUDD 


j k Tjjfc 


NB: One of t 


le following RV 


. . IN blocks must be left out: 


(which one up to user and RGE code) 


RVKAPPAIN 


RVKAPPA 




RVDIN 


RVD 


I Di 


RVSNVEVIN 


RVSNVEV 


i Vi 


RVM2LH1IN 


RVM2LH1 





Table 1: Summary of R-parity violating SLHA2 data blocks. All parameters are given in 
the Super-CKM/PMNS basis. Only 3 out of the last 4 blocks are independent. Which 
block to leave out of the input is in principle up to the user, with the caveat that a given 
spectrum calculator may not accept all combinations. See text for a precise definition of 
the format. 

where the 7 (2-component) generalised neutrinos = i^i, ■■■,^7)'^ are defined strictly 
mass-ordered, i.e., with the i^*^2"'^,3^'^ lightest corresponding to the mass entries for the 
PDG codes 12, 14, and 16, and the four heaviest to the PDG codes 1000022, 1000023, 
1000025, and 1000035 (see also appendix 

Note! although these codes are normally associated with names that imply a specific 
flavour content, such as code 12 being and so forth, it would be exceedingly complicated 
to maintain such a correspondence in the context of completely general mixing, hence we 
do not make any such association here. The flavour content of each state, i.e., of each PDG 
number, is in general only defined by its corresponding entries in the mixing matrix RVNMIX. 
Note, however, that the fiavour basis is ordered so as to reproduce the usual associations in 
the trivial case (modulo the unknown fiavour composition of the neutrino mass eigenstates). 

In the limit of CP conservation, the default convention is that be a real matrix and 
one or more of the mass eigenstates may have an apparent negative mass. The minus sign 
may be removed by phase transformations on x° = z/j as explained in SLHAl [1]. 

Charginos and charged leptons may also mix in the case of L-violation. In a similar 
spirit to the neutralino mixing, we defin^ 

C'^r = -r^M^J-' + h.c. , (47) 

in the basis of 2-component spinors = (e^, fii, tl, —iw~, h^)'^, = {cr, Jir, tr, —iw~^, h2)'^, 
where = (w^ =F w'^)/^/2. Note that in the limit of no RPV the lepton fields are mass 

''Note that the absence of a factor 1/2 on the r.h.s. of eq. (|T7)) corrects and supersedes the pubHshed 
version of this paper. 
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eigenstates. 

We define the unitary 5x5 charged fermion mixing matrices U, V, blocks RVUMIX , 
RVVMIX, such that: 

- tp-^Mj,+^+ = - ip-^U^ U*M^+V^ VtP+ . (48) 



diag(m 



The generahsed charged leptons x~ = (ei, 62, 63, 64, 65) are four-component Dirac fermions, 
and the left-handed and right-handed parts of Cj are the two-component fermions and 
Rf, respectively. They are defined as strictly mass ordered, i.e., with the 3 hghtest states 
corresponding to the PDG codes 11, 13, and 15, and the two heaviest to the codes 1000024, 
1000037. As for neutralino mixing, the flavour content of each state is in no way implied by 
its PDG number, but is only defined by its entries in RVUMIX and RVVMIX. Note, however, 
that the flavour basis is ordered so as to reproduce the usual associations in the trivial 
case. For historical reasons, codes 11, 13, and 15 pertain to the negatively charged field 
while codes 1000024 and 1000037 pertain to the opposite charge. The components of x~ 
in "PDG notation" would thus be (11,13,15,-1000024,-1000037). In the limit of CP 
conservation, U and V are chosen to be real by default. 

R-parity violation via lepton number violation implies that the sneutrinos can mix 
with the Higgs bosons. In the limit of CP conservation the CP-even (-odd) Higgs bosons 
mix with real (imaginary) parts of the sneutrinos. We write the neutral scalars as (j)^ = 
\/2Re |(i?{', H2, 1'e, i^fi, t^r)"^}, with the mass term 

C = -\<P'^Mlo<P\ (49) 

where TW^o is a 5 x 5 symmetric mass matrix. We define the orthogonal 5x5 mixing matrix 
N (block RVHMIX) by 

- 0°^Al5o0° = -£^N7WjoK^K/ , (50) 

diag(m|o) *° 

where = {h^, /i^, h'^, h^, h2^) are the neutral scalar mass eigenstates in strictly increasing 
mass order (that is, we use the label h for any neutral scalar mass eigenstate, regardless 
of whether it is more "Higgs-like" or "sneutrino-like" ) . The states are numbered sequen- 
tially by the PDG codes (25 , 35 , 1000012 , 1000014 , 1000016) , regardless of flavour content. 
The same convention will be followed below for the neutral pseudoscalars and the charged 
scalars. 

We write the neutral pseudo-scalars as 0° = y/2lm ^{H^, H^-, Ve-, ^ii-, ^r)"^}; with the mass 
term 

^ = -\^'^M%-^' , (51) 

where A1|o is a 5 x 5 symmetric mass matrix. We define the 4x5 mixing matrix K (block 
RVAMIX) by 

- 0°^Al|o0° = -l^^M%^^^ , (52) 

diag(m|o) *° 
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where = {A^, A2, A3, A^) are the pseudoscalar mass eigenstates, again in strictly increas- 
ing mass order. The states are numbered sequentially by the PDG codes (36,1000017, 
1000018,1000019), regardless of flavour composition. The Goldstone boson G° (the "5th 
component") has been explicitly left out and the 4 rows of N form a set of orthonormal 
vectors. 

If the blocks RVHMIX, RVAMIX are present, they supersede the SLHAl ALPHA vari- 
able/block. 

The charged sleptons and charged Higgs bosons also mix in the 8x8 mass squared 
matrix Ml± by a 7 x 8 matrix C (block RVLMIX): 



C = - {H^\ Ht, el . e* )Ct CMl^C^ C 



diag(>!2 ) 



/ Hi \ 

Hr 

V J 



(53) 



where ij,k,l e {1,2,3}, a, /3 E {1,...,6} and = 
these states are numbered sequentially by the PDG codes (37, 1000011 , 1000013, 1000015, 
2000011,2000013,2000015), regardless of flavour composition. The Goldstone boson G+ 
(the "8th component") has been explicitly left out and the 7 rows of C form a set of 
orthonormal vectors. 

R-parity violation may also generate contributions to down-squark mixing via additional 
left-right mixing terms. 



(54) 



where Vk are the sneutrino vevs. However, this only mixes the six down-type squarks 
amongst themselves and so is identical to the effects of flavour mixing. This is covered in 
section WAl (along with other forms of flavour mixing). 



4.3 CP Violation 

When adding CP violation to the MSSM model parameters and mixing matrices (for a 
recent review see, e.g., the CPNSH report [44]), the SLHAl blocks are understood to 
contain the real parts of the relevant parameters. The imaginary parts should be provided 
with exactly the same format, in a separate block of the same name but prefaced by IM. The 
defaults for all imaginary parameters will be zero. Thus, for example, BLOCK IMAU, IMAD, 
IMAE, Q= ... would describe the imaginary parts of the trilinear soft SUSY-breaking 
scalar couplings. For input, BLOCK IMEXTPAR may be used to provide the relevant imaginary 
parts of soft SUSY-breaking inputs. In cases where the deflnitions of the current paper 
supersede the SLHAl input and output blocks, completely equivalent statements apply. 

One special case is the fi parameter. When the real part of /i is given in EXTPAR 23, the 
imaginary part should be given in IMEXTPAR 23, as above. However, when \fi\ is determined 
by the conditions for electroweak symmetry breaking, only the phase (f^ is taken as an input 
parameter. In this case, SLHA2 generalises the entry MINPAR 4 to contain the cosine of 
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the phase (as opposed to just sign(/i) in SLHAl), and we further introduce a new block 
IMMINPAR whose entry 4 gives the sine of the phase, that is: 



BLOCK MINPAR 



4 : CP conserved: sign(/i). 
CP violated: cosip^ = 1 



Re{/i}/|/i 



BLOCK IMMINPAR 



4 : CP conserved: n/a. 
CP violated: siny?^ 



Im{/i}/|/i 



Note that cos ip^ coincides with sign(/i) in the CP-conserving case. 

When CP symmetry is broken, quantum corrections cause mixing between the CP- 
even and CP-odd Higgs states. Writing the neutral scalar interaction eigenstates as 0° = 
V2{Re {Hf} , Re {H^} , Im {H^} , Im {H^})^ we define the 3 x 4 mixing matrix S (blocks 
CVHMIX and IMCVHMIX) by 



where $° = {h^,h2,h'^)^ are the mass eigenstates; these states are numbered sequentially 
by the PDG codes (25,35,36), regardless of flavour composition. That is, even though 
the PDG reserves code 36 for the CP-odd state, we do not maintain such a labelling here, 
nor one that reduces to it. This means one does have to exercise some caution when taking 
the CP conserving limit. 

The matrix S thus gives the decomposition of the three physical mass eigenstates in 
terms of the four interaction eigenstates, all in one go, with the Goldstone boson G° ex- 
plicitly projected out and the 3 rows of S forming a set of orthonormal vectors. 

For comparison, in the literature, the projecting-out of the Goldstone boson is often done 
as a separate step, by first performing a rotation by the angle f3. (This is, for instance, the 
prescription followed by CPsuperH [13]). In such an approach, our matrix S would be 
decomposed as: 



where Osxs gives the decomposition of the three physical mass eigenstates in terms of the 
intermediate basis 0° = {^/2Re {H^} , y/2Re {H^} , A^^^J'^ , with A^^^ denoting the tree- 
level MSSM non-Goldstone pseudoscalar mass eigenstate. Note that a simple rotation by 



,0 




(55) 




(56) 
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(3 suffices to translate between the two conventions, so whichever is the more practical can 
easily be used. 

A second alternative convention,e.g. adopted by FeynHiggs [11,38], is to also rotate 
the CP-even states by the angle a as part of the first step. In this case, our matrix S would 
be decomposed as: 



3x3 



/ — sin a cos a ^ \ 

cos a sin a 

— sin /3 cos [3 

V COS/3 sin/3 y 



(57) 



with a defined as the mixing angle in the CP-even Higgs sector at tree-level and 7^.3x3 giving 
the decomposition of the three physical mass eigenstates in terms of the intermediate basis 
$° = {hP , y4°)^gp, that is in terms of the the tree-level mass eigenstates. In order to 
translate between S and 7^.3x3, the tree-level angle a would thus also be needed. This 
should be given in the SLHAl output BLOCK ALPHA: 



BLOCK ALPHA 

CP conserved: a; precise definition up to spectrum calculator, see SLHAl. 
CP violated: atrcc- Must be accompanied by the matrix 5, as described 
above, in the blocks CVHMIX and IMCVHMIX. 

For the neutralino and chargino mixing matrices, the default convention in SLHAl 
(and hence for the CP conserving case) is that they be real matrices. One or more mass 
eigenvalues may then have an apparent negative sign, which can be removed by a phase 
transformation on Xi as explained in SLHAl [1]. When going to CPV, the reason for 
introducing the negative-mass convention in the first place, namely maintaining the mixing 
matrices strictly real, disappears. We therefore here take all masses real and positive, with 
A^, t/, and V complex. This does lead to a nominal dissimilarity with SLHAl in the limit 
of vanishing CP violation, but we note that the explicit CPV switch in MODSEL can be used 
to decide unambiguously which convention to follow. 



5 The Next-to-Minimal Supersymmetric SM 

The first question to be addressed in defining universal conventions for the next-to-minimal 
supersymmetric standard model is just what field content and which couplings this name 
should apply to. The field content is already fairly well agreed upon; we shall here define 
the next-to-minimal having exactly the field content of the MSSM with the addi- 

tion of one gauge-singlet chiral superfield. As to couplings and parameterisations, several 
definitions exist in the literature (for a recent review see, e.g., the CPNSH report [45]). 
Rather than adopting a particular one, or treating each special case separately, below we 
choose instead to work at the most general level. Any particular special case can then be 
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obtained by setting different combinations of couplings to zero. For the time being, how- 
ever, we do speciahse to the SLHAl-hke case without CP violation, R-parity violation, or 
flavour violation. Below, we shall use the acronym NMSSM for this class of models, but we 
emphasise that we understand it to relate to field content only, and not to the presence or 
absence of specific couplings. 

5.1 Conventions 

We write the most general CP conserving NMSSM superpotential as (extending the notation 
of SLHAl): 

Wnmssm = Wmssm - eabXSH'^Hl + ^kS"" + ^/i'^^ + ^pS , (58) 

where Wmssm is the MSSM superpotential, eq. ([2]). A non-zero A in combination with a 
VEV {S) of the singlet generates a contribution to the effective /i term /iefi = A {S) + /i, 
where the MSSM /i term is normally assumed to be zero in NMSSM constructions, yielding 
A^eff = A {S). The sign of the A term in eq. fl58l) coincides with the one in [16,37] where 
the Higgs doublet superfields appear in opposite order. The remaining terms represent a 
general cubic potential for the singlet; k is dimensionless, has dimension of mas^, and 
has dimension of mass squared. The soft SUSY-breaking terms relevant to the NMSSM 

are 

l^soft = V2,MSSM + V:,,MSSM + ml\S\^ + {-ea,\A^SHlHl + \A,S'' + ]^^^^^ , 

(59) 

where V^,mssm are the MSSM soft terms defined in eqs. ([3]) and (jl]), and we have introduced 
the notation = B' ^' . 

At tree level, there are thus 15 parameters (in addition to niz which fixes the sum 
of the squared Higgs VEVs) that are relevant for the Higgs sector of the R-parity and 
CP-conserving NMSSM: 

tan/3, /i, m]j^, mjj^, mj, A, k, Ax, A^, /i', m'|, ^p, ^5, A (5*) , m| . (60) 

The minimisation of the effective potential imposes 3 conditions on these parameters, such 
that only 12 of them can be considered independent. We leave it up to each spectrum 
calculator to decide on which combinations to accept. For the purpose of this accord, 
we note only that to specify a general model exactly 12 parameters from eq. flBUl) should 
be provided in the input, including explicit zeroes for parameters desired "switched off". 
However, since fi = m\ = fi' = m'J = ^i;' = ^5 = in the majority of phenomenological 
constructions, for convenience we also allow for a six-parameter specification in terms of 
the reduced parameter list: 

tan/3, mjj^, m^^. A, k, Ax, A^, X {S) , m| . (61) 

^Note that the factors 1/2 in front of the 11' and m'g terms m eqs. ([55)) and ([55)) . respectively, correct 
and supersede the pubhshed version of this paper. 
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To summarise, in addition to m^, the input to the accord should contain either 12 
parameters from the hst given in eq. fl6Up . including zeroes for parameters not present in 
the desired model, or it should contain 6 parameters from the list in eq. fl61l) . in which case 
the remaining 6 "non-standard" parameters, m|, /i', m'J, ^p, and ^p, will be assumed 
to be zero; in both cases the 3 unspecified parameters (as, e.g., m^^, m'jj^, and m|) are 
assumed to be determined by the minimisation of the effective potential. 

5.2 Input/Output Blocks 

Firstly, as described above in section [3], BLOCK MODSEL should contain the switch 3 with 
value 1, corresponding to the choice of the NMSSM particle content. 

Secondly, for the parameters that are also present in the MSSM, we re-use the corre- 
sponding SLHAl entries. That is, mz should be given in SMINPUTS entry 4 and ■mfj_^,mfj^ 
can be given in the EXTPAR entries 21 and 22. tan /3 should either be given in MINPAR entry 3 
(default) or EXTPAR entry 25 (user-defined input scale), as in SLHAl. If /i should be desired 
non-zero, it can be given in EXTPAR entry 23. The corresponding soft parameter m| can be 
given in EXTPAR entry 24, in the form m|/(cos/3 sin/9), see [1]. The notation m\ that was 
used for that parameter in the SLHAl is no longer relevant in the NMSSM context, but by 
keeping the definition in terms of mg and cos (3 sin /3 unchanged, we maintain an economical 
and straightforward correspondence between the two cases. 

Further, new entries in BLOCK EXTPAR have been defined for the NMSSM specific input 
parameters, as follows. As in the SLHAl, these parameters are all given at the common 
scale Minput) which can either be left up to the spectrum calculator or given explicitly using 
EXTPAR or QEXTPAR (see section 

BLOCK EXTPAR 



Input parameters specific to the NMSSM (in addition to the entries defined in section |2]) 

61 : A (Minput). Superpotential trilinear Higgs SH2H1 coupling. 

62 : K(Minput)- Superpotential cubic S coupling. 

63 : y4A(Minput)- Soft trilinear Higgs SH2H1 coupling. 

64 : (Minput)- Soft cubic S coupling. 

65 : A {S) (Minput)- Vacuum expectation value of the singlet (scaled by A). 

66 : ^^(Mnput)- Superpotential linear S coupling. 

67 : ^^(Minput)- Soft linear 5* coupling. 

68 : /i'( Minput)- Superpotential quadratic S coupling. 

69 : m'J(Minput)- Soft quadratic 5* coupling (sometimes denoted fi'B'). 

70 : m|(Minput)- Soft singlet mass squared. 
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Important note: only 12 of the parameters listed in eq. fl60l) should be given as input at 
any one time (including explicit zeroes for parameters desired "switched off"), the remaining 
ones being determined by the minimisation of the effective potential. Which combinations 
to accept is left up to the individual spectrum calculator programs. Alternatively, for 
minimal models, 6 parameters of those listed in eq. fl6T|) should be given. 

For non-zero values, signs can be either positive or negative. As noted above, the 
meaning of the already existing entries EXTPAR 23 and 24 (the MSSM /i parameter and 
corresponding soft term) are maintained, which allows, in principle, for non-zero values for 
both /i and (S). The reason for choosing A {S) rather than (S) as input parameter 65 is 
that it allows more easily to recover the MSSM limit A, k — ^ 0, (S) oo with A (S) fixed. 

In the spectrum output, running NMSSM parameters corresponding to the EXTPAR 
entries above can be given in the block NMSSMRUN Q= . . . : 



BLOCK NMSSMRUN Q= . . . 

Output parameters specific to the NMSSM, given in the DR scheme, at the scale Q. As in 
the SLHAl, several of these blocks may be given simultaneously in the output, each then 
corresponding to a specific scale. See corresponding entries in EXTPAR above for definitions. 



1 


: A(g)i^i^. 


2 


: k{Q)^. 


3 


: AxiQ)^. 


4 


: A^iQ)^. 


5 


: A {S) {Qf 


6 


: UQ)^. 


7 




8 


: i^'iQ)^. 


9 


: m'|(Q)DR. 


10 


: m|(g)°^. 



5.3 Particle Mixing 

In the CP-conserving NMSSM, the CP-even interaction eigenstates are 
We define the orthogonal 3x3 mixing matrix 5" (block NMHMIX) by 

- 0°^Al5o0° = - 0^5A<5o5^ 5/ , (62) 
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where $° = {h^, /13) are the mass eigenstates ordered in mass. These states are numbered 
sequentially by the PDG codes (25,35,45). The format of BLOCK NMHMIX is the same as 
for the mixing matrices in SLHAl. 

In the MSSM limit (A, k ^ 0, and parameters such that ~ Re {5}) the elements of 
the first 2x2 sub-matrix of Sij are related to the MSSM angle a as 

5*11 ~ — sina , 5*21 ~ cos a , 

5*12 ~ cos a , 5*22 ~ sin a . 

In the CP-odd sector the interaction eigenstates are 

00 = ^\m[{HlHlSf]. We 

define the 2x3 mixing matrix P (block NMAMIX) by 

- 0°^A^^o0° = -l^PMloP^P^ , (63) 

diag(m|o) 

where '^^ = {A^, A2) are the mass eigenstates ordered in mass. These states are numbered 
sequentially by the PDG codes (36,46). The Goldstone boson (the "3rd component") 
has been explicitly left out and the 2 rows of P form a set of orthonormal vectors. An 
updated version NMSSMTools [37] will follow these conventions. 

If NMHMIX, NMAMIX blocks are present, they supersede the SLHAl ALPHA variable/block. 

The neutralino sector of the NMSSM requires a change in the definition of the 4x4 
neutralino mixing matrix to a 5 x 5 matrix. The Lagrangian contains the (symmetric) 
neutralino mass matrix as 

= -Iv^o^A^^oV^o + h.c. , (64) 

in the basis of 2-component spinors ^/^^ = (— —iw^, hi, h2, s)'^. We define the unitary 
5x5 neutralino mixing matrix A^ (block NMNMIX), such that: 

- l^'^M^o^' = ~^;^^N*M^oN^N^ , (65) 

X""^ diag(m^o) 

where the 5 (2-component) neutralinos Xi are defined such that the absolute value of their 
masses increase with i. As in SLHAl, our convention is that A^ be a real matrix. One or 
more mass eigenvalues may then have an apparent negative sign, which can be removed 
by a phase transformation on Xi- The states are numbered sequentially by the PDG codes 
( 1000022 , 1000023 , 1000025 , 1000035 , 1000045) . 



6 Conclusion and Outlook 

At the time of writing of the SLHAl, a large number of computer codes already existed 
which used MSSM spectrum and coupling information in one form or another. This had 
several advantages: there was a high motivation from program authors to produce and 
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implement the accord accurately and quickly, and perhaps more importantly, the SLHAl 
was tested "in anger" in diverse situations as it was being written. 

We find ourselves in a slightly different situation in terms of the SLHA2. There are 
currently few programs that utilise information in any of the NMSSM or CP-violating, 
R-parity violating, or non-trivial flavour violating MSSM scenarios. Thus we do not have 
the beneflt of comprehensive simultaneous testing of the proposed accord and the strong 
motivation that was present for implementation and writing of the original one. What 
we do have are the lessons learned in connection with the SLHAl itself, and also several 
almost-flnished codes which are now awaiting the flnalisation of SLHA2 in order to publish 
their first official releases. Concrete tests involving several of these were thus possible in 
connection with this writeup. 

We have adhered to the principle of backward compatibility wherever feasible. We 
therefore expect that the conventions and agreements reached within this paper constitute 
a practical solution that will prove useful for SUSY particle phenomenology in the future. 
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A PDG Codes and Extensions 



The existing PDG nomenclature for (s)particle names is based on the hmit of the MSSM 
in which CP, R-parity, and flavour are conserved. Several of the mass eigenstates are 
therefore labeled to indicate definite R, CP, and/or flavour quantum numbers. When the 
corresponding symmetries are broken, such a labeling becomes misleading. Throughout 
this paper we have adopted the convention of assigning a common label to all states which 
carry identical conserved quantum numbers in the given model. We then re-use the existing 
PDG codes for those states, arranged in strictly increasing mass order. 

This implies that, while the PDG numbers remain unaltered, their labels change, de- 
pending on which scenario is considered. The PDG codes and labels are discussed in detail 
in the individual sections on flavour violation, R-parity violation, CP violation, and the 
NMSSM. In the tables below, we summarise the PDG numbers and suggested labels rele- 
vant to each distinct scenario, for squarks (Tab. [2]), charged colour-singlet fermions (Tab. [3]), 
neutral colour- singlet fermions (Tab.Hj), charged colour-singlet scalars (Tab. [5]), and neutral 
colour-singlet scalars (Tab. [6l), respectively. Note that these extensions are not officially 
endorsed by the PDG at this time. Codes for other particles can be found in [52, chp. 33]. 



Scalar Quarks 



FLV 


No 


Yes 


RPV 


No 


No 


CPV 


No 


No 


1000001 


dL 


di 


1000002 






1000003 


sl 


d2 


1000004 


cl 




1000005 


h 


d-i 


1000006 


ii 




2000001 


dR 


d4 


2000002 


Ur 




2000003 


sr 




2000004 


Cr 


Up 


2000005 


h 


4 


2000006 




Me 




di 
d2 
ds 
di 

de 



dL 

Ul 

Sl 

Cl 

bi 

ii 

dR 

Ur 

Sr 

Cr 



di 

Ul 

d2 

U2 

d3 
di 

Ui 
Up 

da 



di 

Ul 

d2 

U2 

d, 

Uji 

di 

Ui 

db 

Up 

de 



dL 

Ul 

Sl 

Cl 

bi 

ii 

dR 

Ur 

Sr 

Cr 

h 



Table 2: Particle codes and corresponding labels for squarks. The labels in the first column 
correspond to the current PDG nomenclature. 
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Charged Leptons and Charginos 



FLV 


No 


RPV 


No 


CPV 


No 


11 


e 


13 




15 




1000024 


xt 


1000037 


xt 



Yes I No I No II Yes II Yes I| >^ 
No ^^^^ No II Yes II No II w 
No I I Yes II No II Yes || ^ 



62 

63 
et 



e 

xt 
xt 




e 

xt 
xt 



e 

xt 
xt 



Table 3: Particle codes and corresponding labels for charged colour-singlet fermions. The 
labels in the first column correspond to the current PDG nomenclature. Note that, for 
historical reasons, codes 11, 13, and 15 pertain to negatively charged fields while codes 
1000024 and 1000037 pertain to the opposite charge. 



Neutrinos and Neutralinos 




12 

14 

16 
1000022 
1000023 
1000025 
1000035 
1000045 



No 
No 



No 



xti 
xl 

xl 



IYes_ I No II No II Yes II Yes I| ^ 
No Fyes II No II Yes II No II 1 
No I No II Yes II No II Yes || ^ 



1^1 

xl 

X'2 

xl 
xl 



X? 

xl 
xl 



Z/7 



1^2 

X? 

xt 
xl 
xl 



X? 

x^ 
x^ 
xl 
x% 



Table 4: Particle codes and corresponding labels for neutral colour-singlet fermions. 
labels in the first column correspond to the current PDG nomenclature. 



The 



31 



Charged Higgs Boson and Charged Scalar Leptons 



FLV 


No 




Yes 


sxr V 


No 




No 


CPV 


No 




No 


37 


H+ 


H+ 


1000011 




et 


1000013 


fit 


et 


1000015 


ft 


et 


2000011 




et 


2000013 




et 


2000015 


ft 


et 



ht 


H+ 


ht 


H+ 


H+ 


ht 




ht 


et 


eh 


ht 


fit 


ht 


et 


P't 


ht 


ft 


ht 


et 


ft 


ht 


et 


ht 


et 


et 


ht 




ht 


e^ 


f^R 


ht 


ft 


ht 


e& 


ft 



Table 5: Particle codes and corresponding labels for charged colour-singlet scalars. The 
labels in the first column correspond to the current PDG nomenclature. 



Neutral Higgs Bosons and Scalar Neutrinos 




25 

35 

36 

45 

46 
1000012 
1000014 
1000016 
1000017 
1000018 
1000019 



No 
No 



No 



h' 
A' 



ML 



AO 



1^1 [^is) 



hi 

hi 



hi 
hi 

hi 
Al 

Al 




h\ 
hi 
A\ 



hi 
hi 

hi 
Al 

Al 
Al 



hi 
hi 
hi 



h\ 
hi 

A\ 
hi 
Al 



EL 



Table 6: Particle codes and corresponding labels for neutral colour-singlet scalars. The la- 
bels in the first column correspond to the current PDG nomenclature. The labels in paren- 
thesis denote the optional separation of sneutrinos into separate scalar and pseudoscalar 
components. 
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